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Finite-temperature linear conductance from the Matsubara Green function 
without analytic continuation to the real axis 
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We illustrate how to calculate the finite-temperature linear-response conductance of quantum 
impurity models from the Matsubara Green function. A continued fraction expansion of the Fermi 
distribution is employed which was recently introduced by Ozaki [Phys. Rev. B 75, 035123 (2007)] 
and converges much faster than the usual Matsubara representation. We give a simplified derivation 
of Ozaki's idea using concepts from condensed matter theory and present results for the rate of 
convergence. In case that the Green function of some model of interest is only known numerically, 
interpolating between Matsubara frequencies is much more stable than carrying out an analytic 
continuation to the real axis. We demonstrate this explicitly by considering an infinite tight-binding 
chain with a single site impurity as an exactly-solvable test system, showing that it is advantageous to 
calculate transport properties directly on the imaginary axis. The formalism is applied to the single 
impurity Anderson model, and the linear conductance at finite temperatures is calculated reliably at 
small to intermediate Coulomb interactions by virtue of the Matsubara functional renormalization 
group. Thus, this quantum many-body method combined with the continued fraction expansion of 
the Fermi function constitutes a promising tool to address more complex quantum dot geometries 
at finite temperatures. 

PACS numbers: 71.27.-|-a, 72.10.-d, 73.21.La 



I. INTRODUCTION 

Computing the finite-temperature T linear conduc- 
tance for widely used model systems like an Anderson 
impurity connected to non-interacting leads involves in- 
tegrals over the impurity single-particle spectral func- 
tion p{uj) multiplied by the derivative of the Fermi 
distributionii At T > 0, however, most quantum many- 
particle methods are set up in imaginary time and yield 
the relevant one-particle propagator Q{iuJn) at the Mat- 
subara frequencies a;„ on the imaginary axis. Thus, an 
analytic continuation is required in order to obtain the 
corresponding spectral functions on the real axis and to 
compute transport properties at finite temperatures. If 
Q{iujn) is only known numerically - which is the generic 
scenario - the analytic continuation can in principle be 
achieved, e.g., by Fade approximation of the Matsubara 
data^^i^ but ingeneral such a procedure is ill-controlled 
(see, e.g., Ref.|j). Alternatively, one can use the analytic 
properties of the Fermi function and resort to the calcu- 
lus of residues to perform the integration directly on the 
imaginary axis. Even though this might involve the prop- 
agator Q (or, aiming at the conductance, its derivative) 
not only at the physical Matsubara frequencies but also 
at different arguments, computing such data by interpo- 
lation is expected to be significantly more stable than 
carrying out an analytic continuation. 

In order to rotate the path of integration from the 
real to the imaginary axis, one can employ the Mat- 
subara expansion for the Fermi distribution f{x) = 
l/[exp(a;) + 1], but this is not useful for practical appli- 
cations in fcrmionic many-body problems because of its 



slow convergence!^ A much more efficient representation 
was recently proposed by Ozakii^ It is based on a contin- 
ued fraction expansion and ultimately yields an infinite 
sum over simple poles which arc all located on the imag- 
inary axis. In its infinite form, this series is equivalent 
to the usual Matsubara expansion. Truncation, however, 
modifies the position of the poles (but they remain on 
the imaginary axis) as well as the corresponding residues, 
and it turns out that the resulting finite series yields an 
excellent approximation to f{x) for not too large values 
of I a; I, even if only a few terms are accounted for (all this 
will be quantified below). 

The motivation for the present work is three-fold. 
First, we give a simplified derivation of Ozaki's idea 
which relics on concepts taken from condensed matter 
theory and present results for the rate of convergence of 
the truncated series (Section |lTl. Secondly, we employ an 
exactly-solvable model - a non-interacting tight-binding 
chain with a single site impurity - and demonstrate that 
interpolation of the Matsubara Green function is much 
more stable than continuation to real frequencies, illus- 
trating that it is advantageous to perform calculations 
on the imaginary axis (Section IIII[) . Thirdly, the formal- 
ism is applied to the single impurity Anderson model, 
and the Coulomb interaction U is tackled by an approx- 
imation scheme which (in thermal equilibrium) is most 
easily set up in Matsubara space - the functional renor- 
malization group (FRG)i^ii"— By comparison with very 
accurate numerical renormalization group (NRG) data, 
we demonstrate that the FRG allows for reliably extract- 
ing the finite-temperature linear conductance at small to 
intermediate values of U (Section IIVI) . This was pre- 



viously not possible due to the instability of the ana- 
lytic continuation)^ Whereas within the NRG framework 
the numerical effort grows exponentially when consid- 
ering general models with several correlated degrees of 
freedom and one is generically restricted to situations of 
high symmetry, the Matsubara FRG can easily be gener- 
alized to multi-dot or multi-level quantum dot geometries 
without specific symmetries with only a mild increase in 
computation timci^ Combined with the continued frac- 
tion expansion of the Fermi function, the FRG thus holds 
the promise for future applications addressing transport 
properties of such systems at T > 0. 
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II. CONTINUED-FRACTION 
REPRESENTATION OF THE FERMI FUNCTION 

A. Derivation 

In this Section we give a simplified derivation of 
Ozaki's idea.— The first step is to express the Fermi dis- 
tribution in the form 



/W 



1 



_ 1 / sinh(2;/2) 
T ^ 2 V ^ coshjxji) 



(1) 



where later on a; = /3(a; — /^), with /3 and fi denoting the 
inverse temperature and the chemical potential, respec- 
tively. In order to obtain a continued fraction expansion 
for tanh(a;/2). Ozaki introduced the (generalized hyper- 
geometric) auxiliary functions 



with 0(0) ~ 1 and 

o-{7i) — o.{o, + 1) . . . (a + n — 1) 

It is straightforward to show that 

/I x^\ 
cosh X = F { —, — , sinh x = 
V2' 4 / 



n> I 



3 a;2 

2' 4 



(2) 



(3) 



^F\-..-r] ■ (4) 



The first identity, e.g., is easily confirmed by establishing 
(2n)!/(4"n!) = (l/2)(„) by induction. Thus, tanh(a;) is 
determined by a ratio of the form F{a,x)/ F{a ~ l,x), 
rendering it reasonable to seek for a recurrence relation. 
To this end, one calculates the difference of the power 
series of Eq. ([2]) for a — 1 and a: 



F{a- l,x)- F{a,x) 



a{a — 1) 



F{a + l,x) 



(5) 



Dividing by F(a,x), taking the inverse and replacing a 
by flo -I- n gives 



r„{ao,x) 



+ (ao+«)(ao-l+n)^"+l("0,a;) 



(6) 



FIG. 1: (Color online) Approximation fhiix) to the Fermi 
function obtained from the continued fraction representation 
of Eq. (HH) truncated at different N = 2M. Already for M = 
3, the agreement with the exact result f{x) — l/[exp(a::) + 1] 
is excellent for not too large arguments |a;| < 10. Likewise, 
the A4 — 6 - curve nicely approximates f{x) up to |a;| < 60. 
For larger arguments, all fM{x) tend to 1/2. 



where we have introduced the ratio 



r„(ao,x) 



F{ao + n, x) 
F{ao - 1 + n,x) 



(7) 



If one starts with n = 1, iteration yields an expansion 
for ri(ao,a;) in terms of a continued fraction. From 
Eqs. (m and ^, the corresponding representation for 
[f{x) — l/2]/x is obtained setting oo = 1/2 and replacing 
X by a;^/16. 

For a practical application it is not necessary to write 
down the continued fraction explicitly. Rather, one can 
directly exploit the fact that the very same recursion re- 
lation of Eq. (jni) appears if one aims at computing the in- 
verse of an infinite tridiagonal matrix C associated with 
a semi-infinite one-dimensional tight-binding chain with 
nearest neighbor hopping and starting at site n = 1. Us- 
ing the well-known Fcshbach projection method,— one 
obtains (for n > 1) 



C7 



V ^'^'^-^V n+l.n+1 

(8) 
where C[n] denotes the part of the matrix C projected 
onto the sites greater and equal to n. Comparison with 
Eq. (|ni) and eventually combining with Eqs. ([1]) and (|3]) 
yields 



f{x) = l-j{l\il + zxB)-'\l) 



(9) 



with the only non- vanishing matrix elements of B (which 
decrease like l/(4n) for large n) given by 



B, 



n,n+l 



s, 



n-fl,n 



2v/(2n-l)(2n+l) 



n> 1 . 
(10) 
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FIG. 2: Convergence of the continued fraction expansion of 
the Fermi function to f^{ — 1 for very large negative argu- 
ments a; as a function of the 'chain length' N = 2M. Inset: 
By rescaling with l/n, the four curves collapse, indicating 
that the truncation A/ has to increase only as ^y\x\ in order 
to obtain the same convergence. 



Thus, one finally faces the ordinary eigenvalue problem 



B\K, 



ba\ba) , 



ba e 



(11) 



All ba come in pairs owing to the fact that the state |6q,) 
with components {n\ba) = (— l)"(n|&„) is an eigenvector 
of B with eigenvalue —ba- Then, Eq. ^ can be recast 
as 



m - 1 



(l|6n 



X ^ — ^ 

4^1 + ixba 

a 

Q>0 



Ro 



i/ba 



i/ba 



(12) 



corresponds to a finite termination of the infinite con- 
tinued fraction (see also Ref. Illl ). In order to keep the 
property that all eigenvalues of i? come in pairs {ba,—ba), 
one always truncates at an even number N = 2M of sites 

(M e N). 

The eigenvalue problem posed by Eqs. pH]) and pT|) 
can be solved numerically using standard routines for all 
values M of interest. It turns outi^iiSi that the lowest 60 
percent of the ba are very close to the Matsubara values 
ln{2n — 1)]~^, and the associated residues are approxi- 
mately one. The other inverse eigenvalues increase much 
faster than 7r(2n— 1), and their residues Ra become large. 
As ba G K, all poles of the function /m {x) obtained from 
restricting the sum in Eq. p^ to AI terms are located 
on the imaginary axis. This is an important property 
for the applications discussed in Section HVl which is not 
shared by another recent partial fraction decomposition 
of the Fermi function J^ 

Already for rather small values of M, fnix) yields an 
excellent approximation for f{x) for not too large values 
of |a;|. This is illustrated in Fig. [T] Even at M = 3, the 
approximation /3(a;) to the Fermi distribution is excellent 
up to |a;| < 10. Likewise, f^ix) coincides with f{x) to 
the drawing accuracy for arguments |a;| < 60. For very 
large values of |x|, all functions /m eventually tend to 
fM = 1/2. 

In Fig. [5] we show the convergence of /a/ to one for very 
large negative values of x. Results are shown for four 
different arguments, increasing in absolute value propor- 
tional to 71^. By rescaling with l/n, the four curves col- 



lapse, indicating that A/ has to increase only as \/\x\ m 
order to obtain the same rate of convergence. In passing, 
we note that it might be possible to analytically analyze 
the truncation error along the lines of Ref. [ill 



with Ra = |(1|&q)P/(46^), and a > symbohzing that 
the sum extends over the positive part of the spectrum 
(i.e., all a for which 6q > 0) only. This expansion of the 
Fermi function should be compared to the well-known 
(/3 = I) Matsubara sum decomposition^ 



m - \ 



E 

n=l 



I 



iTr(2n — I) x + iTT(2n — I) 



(13) 

which is slowly converging with w.r.t. n. If the infinite 
matrix B is not replaced by a finite one, Eqs. (jl2p and 
()13p should coincide, implying that Ra ~ I and the pos- 
itive eigenvalues given by ba = l/[7r(2n — 1)]. The fact 
that the latter accumulate at zero is obvious from the 
large n decay of the matrix elements B^n+i- 



B. Convergence 

Ozaki pointed out the enormous acceleration of conver- 
gence one obtains for not too large arguments |a;| by re- 
placing the infinite chain by one of finite length N which 



III. STABILITY OF INTERPOLATION AND 
CONTINUATION OF NUMERICAL DATA 

The linear-response finite-temperature conductance of 
the single impurity Anderson modeU^ is determined 
by an energy integral of the derivative of f{(3io) = 
l/[exp(/3a;) -I- 1] multiplied by the local single-particle 
spectral function p{uj)r^ 



G = -r / p{uj) d^f{l3uj)duj = 2rT ^ i?„ Im 



dQ{iuja) 

(14) 
where we have chose units of ?i = e = 1 (with e being the 
elementary charge) and defined the modified Matsubara 
frequencies uia = T/ba- The prefactor F is a measure for 
the level-lead couplings (see below). The second equality 
of Eq. ^^ follows from rotating to the imaginary axis 
using Eq. P^ - it will be derived explicitly in Section 
IIVI If Q{iuJn) is only known numerically at the phys- 
ical Matsubara frequencies aj„; both p(w) and Q'iiCja) 
need to be determined, e.g., from computing the Fade 
approximation^i^ to Q{iix)n)- Whereas it seems reasonable 
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FIG. 3: (Color online) Derivative of the Matsubara Green 
function associated with an infinite tight-binding chain with 
nearest neighbor hopping t and a single site impurity of 
strength e. Solid lines show the exact result of Eq. (|18|l , sym- 
bols data obtained from computing a Pade approximation to 
Q{iu}„) at different temperatures (and thus different physical 
Matsubara frequencies ujn = (2n -I- l)7rT).— ^i^ The curves at 
e/t = 0.1 and e/t = 4 were both multiplied by a factor of two 
for clarity. 



to expect interpolation (required for extracting G'{iuJa)) 
to be rather stable, an analytic continuation (necessary 
to calculate p(w)) is in general an ill-controlled procedure 
(see, e.g., Ref. |j). In this Section, we address this issue 
explicitly by considering an infinite tight-binding chain 
with a single site impurity as an exactly-solvable test 
system. The latter constitutes a special realization of the 
single impurity Anderson model in absence of Coulomb 
interactions and is associated with a rich local density of 
states (featuring both a continuum as well as one bound 
state). 

The infinite tight-binding geometry is governed by the 
Hamiltonian 



H 



-^E 



H.c, 



eckc. 



o'-o ' 



(15) 



where f > is the hopping amplitude between near- 
est neighbors, and e denotes the strength of an impu- 
rity located at site i = 0. We assume the system to be 
in thermal equilibrium associated with temperature T 
and zero chemical potential. Using standard projection 
techniqueSfiS the local impurity Green function can be 
calculated as 



g{iuj) 



1 



aj>0 



1 



— e — 2t'^gQ{iio) 



iV^ 



(16) 



where go{iuj) is the local propagator of an isolated semi- 
infinite chain: 



5o(J'^) 



1 

2i2 



iuj — isgn(aj) V 4^2 + oj"^ 



(17) 



The imaginary part of the derivative of G (which deter- 
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FIG. 4: (Color online) Tight-binding chain, but now compar- 
ing the exact spectral function of Eq. pO|l with the one ob- 
tained from Pade approximation. The position and spectral 
weight of the bound state is indicated by an arrow. Upper and 
lower panel show results computed with different discretiza- 
tions of the imaginary frequency axis.— 



mines the linear conductance; see Eq. (|14|)) reads 

dGiiuj) _ W (4t2 + t^2 _ g2-) 



Im ■ 



dui 



V4t2 + ^2 (4^2 +UJ^+ e2) 



2 ' 



(18) 



and the retarded Green function can be calculated 
straightforwardly from Eq. (fT6)) by replacing iuj — ?> ui + irj: 



G'^'icj) = <^ '^^- 



1 



sgn.(cj) y/co'-^ —At'^—e-^ir} 



|w| < 2t 
\uj\ > 2t 



(19) 



Thus, the local density of states p{uj) = — ImG"'°*(a;)/7r 
is given by 



p(w) =- 



1 V4t2 - a;2 



7r4i2 



e 

V4i2 -f^ £2 



;Q{\Lu\-2t) 



S(uj- sgn(e)V4i2 + e2) . 



(20) 



It features a continuum as well as a single pole. 

We can now test the stability of both the interpolation 
on the imaginary axis and continuation to real frequen- 
cies of the Matsubara Green function. To this end, we 
compute the Pade approximation (as outlined in Refs.y- 
Q) to G{i(jJn) and from that both Q'{iLJn) as well as p{uj) 
and compare with the exact results of Eqs. (fT8|) and (|20|) . 
respectively. The Fade approximation is calculated nu- 
merically for different meshes w„ = (2n + 1)ttT of Mat- 
subara frequencies determined by the temperature T (as 
well as additional discretization parameters'^). For our 
non-interacting tight-binding chain, both the (exact) lo- 
cal density of states and the derivative Q'{iLu„) do not 
depend on T. 

The results for Q'{iujn) and p{lo) obtained from the 
Fade approximation to the propagator of Eq. (|16p at dif- 
ferent temperatures are shown in Figs. [3] and IH respec- 
tively. It turns out that the former agrees well with the 



exact curve and is in particular independent of T. In con- 
trast, the Pade data for p(w) depends strongly on tem- 
perature (as well as on other discretization parametersi^) 
and is plagued by severe artifacts such as p{ijj) < 0. As 
expected, analytic continuation of a numerically known 
Matsubara Green function is much more unstable than 
interpolating on the imaginary axis. This demonstrates 
the advantage of evaluating physical quantities (e.g., the 
linear conductance which is determined either by G'{iijJa) 
or p{u!)) directly in Matsubara frequency space. 



IV. FUNCTIONAL RG FOR THE SINGLE 
IMPURITY ANDERSON MODEL 

In this Section, wc apply the formalism introduced 
above to compute the linear conductance G of the sin- 
gle impurity Anderson model (SIAM) at finite tempera- 
tures. To this end, wc show how to express G in terms of 
the Matsubara Green function using the continued frac- 
tion expansion of the Fermi function and demonstrate 
that in the exactly-solveable non-interacting (test) case, 
only a few poles (M sa 10) need to be accounted for in 
order to obtain the conductance in agreement with the 
one computed on the real axis. For finite (and not too 
large) Coulomb repulsions, the functional renormaliza- 
tion group - which in equilibrium is most easily imple- 
mented in Matsubara frequency space ~ allows for calcu- 
lating G in good agreement with NRG data. This was 
previously not possible due to the instability of the ana- 
lytic continuation^^ 



A. Model Hamiltonian and linear conductance 
from the imaginary axis 

The Anderson model is governed by the Hamiltoniar>i^ 
H = eY^dU. + U {d\d^ - i) [d\d^ - i) 

+ H ^f" y'lkaCLka + (^RkaCBka) (21) 

k(T 



where Cska and d„ denote fermionic annihilation opera- 
tors associated with (left and right) baths as well as a sin- 
gle impurity, respectively. The latter features an onsite 
energy e and a Coulomb repulsion U between electrons of 
different spin directions a. The baths are assumed to be 
in the wide-band limit leading to an energy-independent 
hybridization F = 27r/9bath^^- 

Given the exact local single-particle spectral function 
p(w), the linear- response conductance of the SIAM at 
temperature T can be calculated aa^^ (in units of ft = 
e = l) 



(22) 
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FIG. 5: (Color online) Linear-response conductance (normal- 
ized to the unitary value Go = l/(27r)) of the non-interacting 
single impurity Anderson model as a function of the level po- 
sition e and at different temperatures T. Solid hues show the 
exact result calculated from the real axis (Eq. (|22|l 'l. symbols 
were obtained from the imaginary axis using the continued 
fraction representation of the Fermi function (Eq. (|23|l '). The 
expansion was truncated with rather small values of M but 
agrees nicely with the exact data nevertheless. 



where f{(3uj) = l/[cxp(/3a;) + 1] denotes the Fermi dis- 
tribution. Employing the expansion of Eq. (J12p . one can 
recast Eq. ([^^ in terms of the derivative of the Matsub- 
ara Green function: 



G = F / f{/3uj)d^piu:)duj 



= F 



Q>0 



dujp{u}) dw 



ttj: 



a>0 






+ c.c. 



xa„ [G'''=*(w)-G^'i^(w)] duj 
2Fr^i?„Im.^^(^"") 



a>0 



dijJo 



(23) 



where we have exploited that the retarded and advanced 
Green functions fjio'/adv ^^^ analytic in the upper and 
lower half of the complex plane, respectively. 



B. Results 



For the non-interacting case, both p(w) and the Mat- 
subara Green function can be computed exactly using 
standard projcctioni^i or equation of motion- techniques. 
They read 



Qiiuj) 



iu — e + i sgn(a;)F 



G = -F / piuj) d^fiMdw 



p{uj) = - 



T 



1 



7r(w-e)2 + r2 ' 



(24) 
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FIG. 6: (Color online) Linear conductance of the single impurity Anderson model at finite Coulomb interactions and different 
temperatures T = 0, T /T — 0.1, T /T = 0.2, and T /V = 0.4 (from top to bottom at e = 0). Solid lines show 'numerically 
exact' data obtained from the NRG framework,— symbols functional RG results computed using the Matsubara FRG scheme 
outlined in Ref. as well as the continued fraction expansion of the Fermi function which underlies Eq. (|23p . Inset to (a): 
Linear conductance as a function of temperature at particle-hole symmetry. Symbols denote FRG data at U /V = 2, the solid 
line a quadratic fit to the latter. The dashed line displays the non-interacting curve. 



and the linear-response eonductanee G as a funetion of 
the level position e and temperature T can be calcu- 
lated (numerically) from Eqs. (P^ and ([25)1 . respectively. 
Comparing both approaches provides another demon- 
stration of the rapid convergence of the imaginary-axis 
expression based on the continued fraction expansion of 
the Fermi function. It turns out (see Fig. [5]) that if the 
series of Eq. (P5)) is truncated with (rather small) values 
M K, 10, the resulting conductance G(e) nicely agrees 
vifith that of Eq. (|22p for arbitrary temperatures from 
T < r to T > r. 

In presence 0/ Coulomb interactions, 'numerically ex- 
act' results for the spectral function of the SIAM can 
be obtained from the NRG frameworkJ^ If one aims at 
describing more complex geometries, however, it is desir- 
able to devise approximate schemes. A recent approach 
is provided by the Matsubara functional renormalization 
groupi which re-formulates a given many-particle prob- 
lem in terms of an infinite set of coupled flow equations 
for single-particle irreducible vertex functions with an in- 
frared cutoff as the flow parameter. Truncation of this 
hierarchy renders the FRG approximate w.r.t. the two- 
particle interaction and can hence a priori be justified 
only in the limit of small U . Despite this fact, applica- 
tion of the most simple truncation scheme ~ which yields 
fiow equations for effective system parameters and can 
thus be regarded as a kind of 'RG enhanced Hartree- 
Fock' approach - to various quantum dot geometries in 
equilibrium turned out to give accurate results for the 
zero-temperature linear conductance even at fairly large 
Coulomb interactionsi^ If one employs a more elaborate 
truncation (which accounts for the frequency dependence 
of the two-particle vertex) , one can in principle compute 
the spectral function of the SIAM in agreement with 
NRG data for small to intermediate values of UA^ How- 
ever, the latter requires an analytic continuation of the 
Matsubara Green functionJ^ which for the SIAM was 



observed to be particularly ill-controlled if both T 7^ 
and e y^ Oi^ For this reason, it was previously not pos- 
sible to address the linear conductance as a function of 
the level position at finite temperatures. The continued 
fraction expansion for the Fermi function now allows for 
calculating G(e, T) directly from the imaginary axis. To 
this end, we extract the FRG approximation^^ to the 
Matsubara Green function of the SIAM using precisely 
the formalism outlined in Ref.yj. Thereafter, we compute 
the derivative G'{ii^a) from interpolating by virtue of a 
Fade approximation - which again turns out to be far 
more stable than continuation to the real axis - and ul- 
timately the linear conductance by Eq. ((23)) . The result 
is shown (and compared to NRG datapi^) in Fig. [S] The 
FRG correctly describes the widening of the Lorentzian 
lineshape of G(e) due to Coulomb correlations (signal- 
ing the development of a Kondo plateau) and for small 
U/T = 2 (intermediate U/T = 4) agrees nicely (decently) 
with the NRG reference. Moreover, the temperature- 
dependence of G{T) at particle-hole symmetry e = is 
purely quadratic (see the inset to Fig. [Hlja)) as expected 
from Fermi-liquid theoryj^ For larger values of U, the 
agreement with the numerical RG data becomes worse 
(which is in line with the results of Ref. [J) . 



V. CONCLUSION 

In this Paper we illustrated how to compute the finite- 
temperature linear-response conductance G of quantum 
impurity models from the Matsubara Green function 
Q{iiOn) using a rapidly-converging continued fraction ex- 
pansion of the Fermi function recently introduced by 
Ozaki. In case that Q{iujn) is only known numerically, 
this formalism allows to circumvent the need to carry out 
an (ill-controlled) analytic continuation to the real axis. 
As an application, we studied the single impurity Ander- 



son model at finite temperatures within the framework 
of the Matsubara functional renormalization group and 
showed that G can be obtained accurately in compari- 
son with numerical RG data for not too large Coulomb 
interactions. The FRG therefore holds the promise for 
future treatments of more complex quantum dot geome- 
tries which cannot be easily addressed within the NRG 
framework. 
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